Degree heterogeneity in spatial networks with total cost constraint 
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Recently, In [Phys. Rev. Lett. 104, 018701 (2010)] the authors studied a spatial network which is 
constructed from a regular lattice by adding long-range edges (shortcuts) with probability Pij ~ ^Tj"', 
where rtj is the Manhattan length of the long-range edges. The total length of the additional 
edges is subject to a cost constraint (^r — C). These networks have fixed optimal exponent a 
for transportation (measured by the average shortest-path length). However, we observe that the 
degree in such spatial networks is homogenously distributed, which is far different from real networks 
such as airline systems. In this paper, we propose a method to introduce degree heterogeneity 
in spatial networks with total cost constraint. Results show that with degree heterogeneity the 
optimal exponent shifts to a smaller value and the average shortest-path length can further decrease. 
Moreover, we consider the synchronization on the spatial networks and related results are discussed. 
Our new model may better reproduce the features of many real transportation systems. 

PACS numbers: 89.75.Hc, 02.50.-r, 05.40.Fb, 89.75.Fb 



Spatial features play a significant role in the trans- 
portation networks Internet 0, mobile phone net- 
works 3, power grids [1], social networks [H and neural 
networks !6|. In the past decade, many systems have been 
modeled by complex networks where nodes and links are 
embedded in space. In these models, nodes are located 
in the plane and the geometric distance between nodes 
is well defined. Then links are constructed according to 
rules based on spatial indices J? -[l2| . Growing spatial net- 
works were also studied [H, [ij] . Generally, these model 
of spatial networks are able to reproduce some properties 
of the real systems such as community structure, scale- 
free connection length distribution and so on. For a detail 
review of the field, see [l^ . 

The most important features of spatial network model 
is that there is a cost associated with the length of links, 
which has dramatic effects on the topology and function 
of these networks. A total cost constraint has been re- 
cently introduced to design the spatial networks jT6l - l2l| . 
The total cost C is defined as the total length of the links, 
C where r is the length of links. In [l^, pairs of 

sites ij in 2-dimensional lattices are randomly chosen to 
be linked with long-range connections with probability 
Pij ~ fij" J where is the Manhattan distance between 
sites i and j. New links are added until the total length 
of the links reach the total cost C. The exponent a con- 
trols the trade-off between the link length and link num- 
ber. A large value of a allows for the formation of many 
short-length links while the small a favors the creation 
of a few long-length links. The authors in [l^ show that 
the optimal exponents for both the average shortest-path 
length and navigation steps are a = 3 in 2-dimension 
and a = 2 in 1-dimension (the 1-dimensional scenario 
has been strictly proved very recently in [l^). The au- 
thors claim that such model reveals the optimized aspect 
of airline networks under the conditions of geographical 
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availability (for customer satisfaction) and cost limita- 
tions (for airline company profit). Similar works have 
been carried out to study dynamics such as traffic con- 
gestion and synchronization on spatial networks under 
total cost constraint, a.nd the optimal link length expo- 
nents arc found [H, 0, [HI . 

However, we observe that the degree distribution in 
these spatial network models with total cost constraint 
is homogenous, which is far different from real cases, in- 
cluding airline systems. Many previous works have re- 
vealed that many geographic-based networks have het- 
erogenous degree distribution [l|. In this paper, we pro- 
pose a method to introduce degree heterogeneity in spa- 
tial networks with total cost constraint. We find that 
degree heterogeneity affects the optimal exponent for the 
link length distribution. Specifically, the exponent shifts 
to a smaller value. Also, the average shortest-path length 
can further decrease due to degree heterogeneity. More- 
over, we study the effect of degree heterogeneity on the 
dynamics taking place on the spatial networks, and we 
observe the same shifting phenomenon of the optimal ex- 
ponent. 

To begin our analysis, we first briefly describe the orig- 
inal spatial network model introduced in [l^ . Nodes are 
located in a d-dimensional regular square lattice, where 
each site i is connected with its 2c? nearest neighbors. 
Pairs of sites ij are randomly chosen to receive long-range 
connections with probability proportional to t",^", where 
rij is the Manhattan distance between sites i and j (i.e., 
the number of connections separating the nodes in the 
underlying regular lattice). Long-range connections are 
added to the system until their total length (cost) ^ rij 
reaches a given value C. The exponent a controls the av- 
erage length of the long-range connections, a smaller a is 
corresponding to fewer but longer connections, and vice 
versa. In this model, the most interesting phenomenon is 
that the average shortest-path length and the navigation 
steps of the spatial networks can be minimized under a 
certain value of the parameter a. Specifically, optimal 
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ar , where a is determined from the normahzation 
condition X^r"^ ^('") ~ ^■ 

3. A node i is chosen randomly, and one of the Nr 
nodes who are at distance r from node i is picked 
randomly, and an edge between them is created if 
there is no edge between them yet. 
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FIG. 1. (Color online) The average shortest-path length {1} 
as a function of 5 in spatial networks [l^ embedded in (a) 
1-dimensional lattices and (b) 2-dimensional lattices, (c) and 
(d) are respectively the degree distribution of these network 
models (N = 1000 for 1-dimension and N = 10000 for 2- 
dimension, S — 2) and the degree distribution of the air- 
line network of USA. Results for the 1-dimensional and 2- 
dimensional models are averaged over 100 and 10 independent 
realizations, respectively. 

values are found to be a* = 3 in 2-dimensional spatial 
networks and a* = 2 in 1-dimensional ones. Accordingly, 
it is concluded that the best transportation condition is 
obtained with an exponent a* = d + 1 where d is the 
dimension of the underlying lattice. 

The probability for two nodes {i and j) to receive a 
long-range connection can be described by the probabil- 
ity density function (PDF) p(r) of long-range connections 
with Manhattan length Actually, p{r) is the summa- 
tion of Pij with rij = r. In d-dimensional space, the 
number of nodes which have Manhattan distance r from 
a given node is proportional to r'^~^ . Therefore, p{r) can 
be expressed as 

p(r) ~ r'^^V-" = r''^"-\ (1) 

which means p{r) ^ r^^ and S = a + 1 — d. Since the 
optimal exponent is a* = d + 1, (5* = (d + 1) + 1 — d = 2 
independently of spaces' dimension. 

Based on this observation, the spatial network in [l^ 
can be equivalently generated by the following procedure: 

1. N nodes are arranged in a d-dimensional lattice 
with periodic boundary condition. Every node is 
connected with its nearest neighbors, which can 
keep every node reachable. In addition, between 
any pair of nodes there is a well defined Manhat- 
tan distance. 

2. A certain distance r (2 < ?' < J'max, where Tmax is 
the largest distance between any nodes in the ini- 
tial network) is generated with probability P{r) — 



4. Update the total cost ^ r. Repeat step 2 and 3 un- 
til the total cost reaches C — cN (In the following, 
we fix c = 10 for convenience). 

We first study the average shortest-path length (l) in 
spatial networks generated by the above procedure (see 
fig. 1(a) and (b)). As expected, the optimal (l) appears 
under the exponent S — 2 in both 1-dimension and 2- 
dimension networks, and this feature is independent of 
the network size. In addition, we observe that the degree 
in such spatial networks is homogenously distributed, as 
shown in fig. 1(c). However, the degree distribution in 
real system is far more heterogeneous [l5| . For example, 
the Chinese airline JE^], the Italian airline [2^ and the 
world wide airline j24| networks are found to have trun- 
cated power-law degree distribution. We show in fig. 1(d) 
that for the airline network of USA (USAir) [l^ the de- 
gree distribution approximately obeys a power-law. To 
capture this feature of real networks, it is necessary to 
introduce degree heterogeneity in spatial network mod- 
els and investigate how the topology and function of the 
networks are affected. Accordingly, we propose a revised 
spatial network generating procedure. Instead of ran- 
domly adding the long connection to the underlying lat- 
tice in step 3, we first calculate a score Sij for each pair 
of node ij with distance r (generated in step 2) and add 
a link between the pair ij with probability proportional 
to its score. Sij is obtained by 

Sij — (kikj) , (2) 

where is a tuneable parameter. When 9 = 0, the gener- 
ated network reduces to the original model in [l^ . When 
< 0, new links will connect the nodes with lower degree, 
which leads to an even more homogeneous degree distri- 
bution. When 9 > 0, new links will preferentially attach 
to the nodes with higher degree, and a heterogenous de- 
gree distribution emerges. Two examples are shown in 
fig. 2(a) and (b). Obviously, by adjusting the parame- 
ter 9, the degree distribution is no longer Poissonian as 
in the original model. When 9 is small (e.g. 9 = 2), the 
spatial networks exhibit a power-law-like degree distribu- 
tion. In this case, the network has several sub hubs which 
are connected to each other by some long distance links. 
When 9 is relatively large (e.g. 9 = 5), a few super hubs 
will emerge and most of the links are connected to such 
nodes. Therefore, the tail of the degree distribution will 
become longer and longer as 9 increases. The degree het- 
erogeneity can be measured by the index H = 
The higher the H, the larger the degree heterogeneity. As 
shown in fig. 2(c), the index H increases with 9, which 
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FIG. 2. (Color online) The degree distribution of the (a) 1- 
dimensional and (b) 2-dimensional modified spatial network 
model under diflerent 6. (c) shows the index H as a. function 
of the parameter 6 and (d) is the link length distributions for 
different 9 in 1-dimensional spatial networks. The parameters 
are set as 5 = 2, = 1000 for 1-dimension and = 10000 for 
2-dimension. Results for the 1-dimensional and 2-dimensional 
models are averaged over 100 and 10 independent realizations, 
respectively. 



confirms that the parameter 6 can indeed control the de- 
gree heterogeneity. In the USAir network, H = 3.46. It 
suggests that the degree heterogeneity of real system is 
far from that of the original network model where 9 = Q 
and H » 1. Furthermore, the modified spatial model 
can perfectly preserve the link length distribution as the 
original model since it doesn't affect the link length gen- 
erating process but only rearranges the location of the 
links. We can see from fig. 2(d) that the link length 
distributions for different 6 is identical as for the case 
= 0. 

Since the average shortest-path length (l) refiects the 
general transportational ability for a network, we con- 
duct simulations for different values of 9 and see how the 
degree heterogeneity affects the optimizing process of (/) 
based on different 6. The simulations are carried out 
in both 1-dimensional and 2-dimensional spatial models 
and the results are reported in fig. 3. When 9 < 0, 
the minimal (l) is always achieved at S* = 2. However, 
when the degree distribution becomes heterogenous (i.e., 
9 > 0), the optimal exponent 6* will shift to a smaller 
value. The results are similar in both 1-dimensional and 
2-dimensional spatial networks (see fig. 3(a) and (b) re- 
spectively) , and imply that a few more longer connections 
are needed to minimize (l) when network degree distri- 
bution is heterogenous. 

For each 9 there is an optimal S* and a minimal (Z*), 
another interesting aspect is hence studying the relation 
between 9 and (/*). As we can see from fig. 3(c) and 
(d) , introducing the degree heterogeneity can indeed fur- 
ther decrease (l*). However, if one keeps enlarging 9 , 
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FIG. 3. (Color online) The dependence of (l) on S under dif- 
ferent 8 in spatial networks generated from (a) 1-dimensional 
lattices with N = 1000 and (b) 2-dimensional lattices with 
N = 10''. The white lines mark the optimal 5* under different 
9. (c) and (d) are the relation between 9 and the minimal (l*) 
in 1-dimensional and 2-dimensional spatial networks respec- 
tively. The results for the 1-dimensional and 2-dimensional 
models are averaged over 100 and 10 independent realizations 
respectively. 



(l*) will start to increase. We can see from fig. 3(c) and 
(d) that optimal 9* for {I*) exist in both 1-dimensional 
and 2-dimensional spatial networks. Interestingly, these 
optimal 9* is independent of the network size but depen- 
dent on the network dimension. Specifically, 9* ~ 2 in 1- 
dimensional spatial networks and 9* = 3 in 2-dimensional 
spatial networks, which leads us to conjecture that the 
optimal value is obtained at 9* = 1 + d. We also con- 
sider the greedy navigation algorithm used in ref. (l6j 
to model the case when only local information of the 
network is known during the transport process. Similar 
optimal exponent shifting phenomenon is observed and 
the minimum navigation steps can further decrease due 
to the degree heterogeneity. 

Another important aspect to consider in spatial net- 
works is how the network function (dynamics) can be 
enhanced. Specifically, refs. HO, HH have revealed that 
there is an optimal link length distribution for the net- 
work synchronizability. Under the framework of mas- 
ter stability analysis, the synchronizability of an undi- 
rected network can be quantified by the eigenvalue ratio 
of the corresponding Laplacian matrix of this network, 
namely R = where Xn and A2 arc respectively 

the largest and the smallest non-zero eigenvalues of the 
Laplacian matrix [26| . Generally, the smaller the value of 
R, the stronger the network synchronizability. In ^0ll2]| . 
starting from 1-dimensional lattice, the authors find the 
optimal length distribution of links (power-law distribu- 
tion with exponent S = 1.5) to add to the network for 
minimizing R. However, these previous network mod- 
els also yields homogeneous degree distribution. Here we 
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FIG. 4. (Color online) (a) The dependence of the optimal 
S* for synchronizability R on 6. (b) The optimal synchroniz- 
ability R* as a function of 9. The spatial networks are built 
on 1-dimensional lattices. The results are averaged over 100 
independent realizations. 



are interested in how the degree heterogeneity affects the 
synchronizability for spatial networks. The results are 
reported in fig. 4. Obviously, the optimal S* for R shifts 
from 6* = 1.5 to a smaller value as 9 increases. However, 
the optimal synchronizability R* cannot be further en- 
hanced by introducing the degree heterogeneity as shown 
in fig. 4(b). It is well known that degree heterogene- 
ity will weaken the synhcornizability while a small av- 
erage shortest-path length is generally favorable for syn- 
chornizability [26|. As we discussed above, introducing 
the degree heterogeneity can further decrease the average 
shortest path. However, the decreased average shortest- 
path length in current case is not sufficient to overcome 
the effect of degree heterogeneity, and hence R* keeps 
increasing with 9. 



To summarize, we have proposed a method to intro- 
duce degree heterogeneity in spatial networks without 
destroying the power-law link length distribution of the 
original model [l^. We have investigated its effect on 
network topology and function. When optimizing the 
average shortest-path length and the network synchro- 
nizability, we find that the best power-law exponent for 
the link length distribution shifts to a smaller value than 
that of the original model. Additionally, the average 
shortest-path length can further decrease due to the de- 
gree heterogeneity while the network synchronizability is 
weakened if the degree distribution is heterogenous. 

From the practical point of view, the original spatial 
network model can explain the power-law connec- 
tion length distribution with given exponent in some real 
systems. By building the relation between the optimal 
exponent and the degree heterogeneity, our model can 
extend such explanation of real transportation networks 
to a wider range of connection length distribution which 
is found to be existing by many empirical studies p^ . 
Furthermore, since we reveal the optimal degree hetero- 
geneity for average shortest path, our model may ex- 
plain why the degree heterogeneity is widely exhibited 
in real systems and only within certain range of values. 
In this sense, our model can not only better understand 
and model real transportation systems, but also help to 
design efficient ones. 
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